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SUMMARY 

Anisotropic high-temperature alloys are used to meet the safety and dur- 
ability requirements of turbine blades for high-pressure turbopumps In reusable 
space propulsion systems. This study assesses the applicability to anisotropic 
components of a simplified Inelastic structural analysis procedure developed 
at the NASA Lewis Research Center. The procedure uses as Input the history of 
the total strain at the critical crack Initiation location computed from elas- 
tic finite-element analyses. Cyclic heat transfer and structural analyses 
were performed for the first stage high-pressure fuel turbopump blade of the 
space shuttle main engine. The blade alloy Is directionally solidified MAR-M 
246 (nickel base). The analyses were based on a typical test stand engine 
cycle. Stress-strain histories for the airfoil critical location were computed 
using both the MARC nonlinear finite-element computer code and the simplified 
procedure. Additional cases were analyzed In which the material yield strength 
was arbitrarily reduced to Increase the plastic strains and, therefore, the 
severity of the problem. Good agreement was shown between the predicted stress 
strain solutions from the two methods. The simplified analysis used about 
0.02 percent (5 percent with the required elastic finite-element analyses) of 
the CPU time used by the nonlinear finite element analysis. 


INTRODUCTION 

Hot gas-path components of advanced aircraft gas turbine engines and 
rocket engines for reusable space propulsion systems operate under extreme gas 
pressure and temperature environments. These operating conditions subject the 
high-pressure stage turbine nozzles and blades to severe thermal transients 
that can result In large Inelastic strains and rapid crack Initiation. To 
attain the safety and durability requirements for these components frequently 
necessitates the use of advanced turbine blading alloys. These alloys exhibit 
mechanical property anisotropy. Assessing or Improving the durability of hot 
section components Is contingent on accurate knowledge of the stress-strain 
history at the critical location for crack Initiation. 

Nonlinear finite-element analysis techniques have become available In 
recent years for calculating Inelastic structural response under cyclic load- 
Ing. Thesennethods~are based~on- classical- Incremental-plastlclty-theory-wlth 
uncoupled creep constitutive models. Many of the nonlinear finite-element 
computer codes such as MARC (ref. 1) have the capability of handling materials 
with anisotropic properties. However, these codes are usually too costly and 
time consuming to use In the early design stages for aerospace applications. 
Costs are further Increased by the geometrical complexity of high-pressure 
turbine blades which require three-dimensional analyses and sometimes substruc- 
turing to obtain accurate solutions. To Improve the design of engine hot path 
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components such as turbine blades, simplified and more economical procedures 
for representing structural response under cyclic loading have been under 
development' (refs. 2 to 4). 

The objective of this study was to evaluate the utility of a simplified 
cyclic structural analysis method In calculating the local stress-strain 
response of an anisotropic turbine blade airfoil at the critical location for 
crack Initiation. The first high-pressure stage fuel turbine blade (HPFTB) In 
the liquid hydrogen turbopump of the space shuttle main engine (SSME) was 
selected for this study. In the past these blades have undergone cracking In 
the blade shank region and at the airfoil leading edge adjacent to the plat- 
form. To achieve the necessary durability, these blades are currently being 
cast using directional solidification. Single crystal alloys are also under 
Investigation for future SSME applications. MARC elastic and elastic-plastic 
finite-element analyses were performed for the blade airfoil. Because of the 
extensive computation time required for the nonlinear finite-element analyses, 
neither the blade platform nor shank regions were modeled. The history of the 
total strain calculated at the critical location from the elastic finite- 
element analysis was used as Input for the simplified procedure. Solutions 
from the simplified Inelastic analyses of these problems for the critical 
airfoil location were compared to those from the MARC nonlinear analyses. 


PROBLEM DESCRIPTION 

The airfoil of the high pressure stage turbine blade of the SSME fuel 
turbopump was analyzed because of Its history of early crack Initiation. This 
blade Is Illustrated In figure 1. The uncooled airfoils have a span length of 
2.2 cm and a span-to-chord width aspect ratio of approximately unity. The 
blades are directionally cast from MAR-M 246+Hf alloy. Temperature-dependent 
properties for this alloy were mainly provided by the Rocketdyne Division of 
Rockwell International Corporation. Material elastic properties are summarized 
In table I. Mean thermal coefficient of expansion data were converted to 
Instantaneous values for MARC Input. Longitudinal stress-strain properties, 
summarized In table II, were used for the elastic-plastic region; transverse 
stress-strain properties were not available at the time of this study. A sin- 
gle crystal alloy Is also being considered for turbine blades In future SSME 
applications. 

Cracking has occurred during service at the airfoil base near the leading 
edge and In the blade root shank area. These cracks were apparently Initiated 
during the first few mission cycles due to the severe thermal transients and 
were propagated by vibratory excitation. Since the primary purpose of this 
study was to compare nonlinear finite-element and simplified analytical methods, 
the blade root and platform were excluded from the analysis to limit the size 
of the problem and, therefore, reduce the computing time. 

The mission used for the analysis Is shown In figure 2 In terms of turbine 
Inlet temperature, gas pressure and RPM. This cycle Is applicable to a factory 
test of the engine; It Is also reasonably representative of a flight mission 
except for the foreshortened steady-state operating time. The major factor 
Inducing fatigue cracking Is the transient thermal stressess caused by the 
sharp Ignition and shutoff transients. 
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Transient and steady-state three dimensional heat transfer analyses have 
been conducted using the MARC code. Film coefficients were obtained from pre- 
liminary Information supplied by Rocketdyne. The gas temperature was assumed 
constant around the airfoil surface for each time step. Colder boundary con- 
ditions were assumed at the airfoil base to simulate the effects of the cooling 
of the blade-to-dlsk attachment region by the liquid hydrogen fuel. 


ANALYTICAL PROCEDURE 

Elastic-plastic analyses have been conducted for the HPFTB airfoil with 
both a simplified analytical procedure developed at the NASA Lewis Research 
Center and with the MARC code. The severity of the problem was progressively 
Increased by analyzing a series of cases In which the material yield strength 
was arbitrarily reduced until plastic strain reversal was obtained In the cycle. 
Separate HARC analyses were conducted for one case using both orthotropic elas- 
tic constants and the Young's modulus and Poisson's coefficient with respect 
to the longitudinal (spanwlse) direction; this was to determine If the longi- 
tudinal properties would give a sufficiently accurate elastic-plastic finite- 
element solution to be used for the simplified analysis. However, the best 
results were obtained with the simplified procedure by the use of effective 
elastic moduli based on MARC elastic finite-element analyses with orthotropic 
material properties. Creep analyses were not conducted because the combination 
of airfoil temperatures and mission dwell times were not severe enough to Induce 
a significant creep problem. Also, there was Inadequate knowledge of the creep 
characteristics of the anisotropic material to perform such analyses even If 
desired. 


Simplified Analysis 

The simplified analytical procedure was developed to economically calcu- 
late the stress-strain history at the critical fatigue location of a structure 
subjected to cyclic thermomechanical loading. This procedure has been exer- 
cised on a wide variety of problems Including multlaxlal loading, nonlsothermal 
conditions, different materials and constitutive models, and dwell times at 
various points In the cycles. Comparisons of the results of the simplified 
analyses with MARC Inelastic solutions for these problems have shown reasonably 
good agreement (refs. 2 to 4). 

The basic assumption Is that the total strain ranges calculated from 
linear elastic and nonlinear Inelastic analyses are approximately equal and, 
therefore, the material cyclic response can be calculated using as Input the 
total strain history obtained from an elastic analysis. This assumption Is 
essentially true for thermally dominated loading. There Is a version of the 

procedure that~uses Neuber-type~cor recti ons-to account-for-straln-redTst-r-Ibu 

tlon due to mechanical loading; however, this version was not utilized for 
this study because of the dominance of the thermal loading during the peak 
strain parts of the cycle. Classical Incremental plasticity methods are used 
to characterize the yield surface by a yield condition to describe yielding 
under multlaxlal stress states and by a hardening model to establish the loca- 
tion of the yield surface during cycling. This procedure can accommodate 
Itself to any yield criterion or hardening model. The only requirements are 
that the elastic Input data, whether calculated or measured, be In a form 
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consistent with the yield criterion and that the appropriate material 
properties be used In conjunction with the hardening model. 

In these analyses, a bilinear kinematic hardening model was used to repre- 
sent the effect of cycling on the yield condition. Since the cyclic stress- 
strain relation Is a function of the plastic strain range, It Is necessary to 
Iterate between the Initially assumed and the calculated maximum plastic 
strains. This Iterative process Is usually accomplished within three Itera- 
tions. However, each Iteration results In some change In the size and shape 
of the cyclic stress-strain loop. These changes, although generally small, 
create some difficulty In directly comparing solutions from the simplified 
procedure against finite-element analysis results because of the differences 
In the stress-strain curves. 

As In most nonlinear computer codes, the von Mlses yield criterion has 
been used In applying stress-strain results from elastic finite-element analy- 
ses of multlaxlal problems as Input for the simplified procedure. To compute 
cyclic hysteresis loops for life prediction purposes, the Input von Mlses 
stresses and strains have to be assigned signs, usually on the basis of the 
signs of the dominant principal stresses and strains. 

The elastic Input data are subdivided Into a sufficient number of Incre- 
ments to define the stress-strain cycle. As will be discussed later, elastic 
finite-element solutions for 6 points In the SSME mission cycle proved adequate 
as the basis to create the total strain history required as Input for the sim- 
plified analysis of the HPFTB airfoil. These points were at the start and end 
of the mission and at the maximum and minimum temperature peaks during the 
preignition and main Ignition phases. A total of 120 stress-strain-temperature- 
time Increments were obtained by Interpolation from the 6 elastic finite- 
element solutions for the critical location. These Increments are analyzed 
sequentially to obtain the cumulative plastic and creep strains and to track 
the yield surface. 

An Iterative procedure Is used to calculate the yield stresses for 
Increments undergoing plastic straining. First, an estimated plastic strain 
Is assumed for calculating an Initial yield stress from the stress-strain 
properties and the simulated hardening model. Then a new plastic strain Is 
calculated as the difference between the total and elastic and creep strain 
components. The yield stress Is then recalculated using the new plastic strain. 
This Iterative procedure Is repeated until the new and previous plastic strains 
agree within a tolerance of 1 percent. Creep computations are performed for 
Increments Involving dwell times using the creep characteristics Incorporated 
In the code. Depending on the nature cf the problem, the creep effects are 
determined on the basis of one of three options to be selected; (1) stress 
relaxation at constant strain, (2) cumulative creep at constant stress, or (3) 
a combination of (1) and (2). 

A FORTRAN IV computer program (ANSYMP) was created to automatically 
Implement the simplified analytical procedure. A detailed description of the 
calculatlonal scheme Is presented In previous papers (refs. 2 to 4) on the 
development of this procedure. 
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MARC Finite-Element Analysis 

A three-dimensional finite-element model of the airfoil (fig. 3) was con- 
structed of eight-node Isoparametric elements. The model consisted of 360 
elements with 576 nodes and 1661 unsuppressed degrees of freedom. The blade 
base and most of the platform were omitted for the MARC nonlinear analysis to 
reduce the computing time and to run the problem In-core on the CRAY computer 
system at Lewis. Boundary conditions were applied to constrain all nodes at 
the base of the model to lie on a platform plane. Additional boundary condi- 
tions were Imposed to prevent rigid body motion In this plane. 

The MARC code has been used extensively at NASA Lewis for Inelastic analy- 
ses of aircraft turbine blades and combustor liners and of space power compo- 
nents. In conducting a cyclic analysis, the loading history Is divided Into a 
series of Incremental load steps which are sequentially analyzed. The plasti- 
city algorithm Is based on a tangent stiffness approach In which the stiffness 
matrix Is reformulated and reassembled for every plastic load Increment. The 
Incremental loads are modified by residual load correction vectors to Insure 
that the solution does not drift from a state of equilibrium. Convergence for 
the Iterative plasticity analysis Is Indicated when the strain energy used In 
assembling the stiffness matrix approximately equals the energy change result- 
ing from the Incremental solution. 

The temperature tolerance controls on the MARC transient heat transfer 
analysis resulted In the automatic subdivision of the mission cycle Into 124 
time Increments. The same Increments were used for the elastic-plastic struc- 
tural analyses. Incremental loading Included centrifugal and gas pressure 
loads and metal temperature distributions as calculated from the heat transfer 
analysis. Approximately one million words of core storage on a CRAY-1S com- 
puter were needed to run the problem. Each cycle of analysis required about 
3 hr of central processor unit (CPU) time on the CRAY system. In terms of 
calendar time, the situation was even more serious because the system was so 
heavily loaded that such a large block of computing time normally was only 
available over weekends. 

The directionality of the elastic material properties causes anisotropic 
constraints. Lekhnltskll (ref. 5) has derived the generalized elastic strain 
equations for an anisotropic body with a transverse plane of Isotropy. Matrix 
Inversion of these equations to solve for the stresses results In the 
relationship 
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where n = E/E' and a = nE'/((l + v)(l - v -2 nv' 2 )). Here E ' , G', and 
v' denote the Young's modulus, shear modulus and Poisson's ratio, respectively, 
for the longitudinal or span direction while E, G, and v denote these con- 
stants with respect to any direction In the transverse plane of Isotropy. 
Rocketdyne supplied values of 0 .143 and 0.391 were used for v' and v, 
respectively. This anisotropic stress-strain law was Incorporated In the MARC 
user subroutine, H00KLW. 

Plastic strain calculations were based on Incremental plasticity theory 
using the von Mlses yield criterion, the normality flow rule and a kinematic 
hardening model. The material elastic-plastic behavior was specified by the 
yield strengths and work hardening properties In the longitudinal direction; 
transverse properties were not available. 


DISCUSSION OF RESULTS 

Calculated metal temperatures at the leading edge at midspan and at the 
crack Initiation site at the base of the airfoil (critical location) are pre- 
sented In figure 4 as a function of elapsed time during the cycle. The assumed 
gas temperature around the airfoil Is also Indicated. Of particular note Is 
that the leading edge temperature at the airfoil base Is cooler than at mldspan 
throughout the cycle. This seems reasonable because of the cooling of the 
blade-to-dlsk attachment region by the liquid hydrogen fuel. The colder air- 
foil base temperatures Induce tensile thermal stresses at the critical leading 
edge location that are additive to the centrifugal stresses. 

The entire discussion of the structural analysis results for the HPFTB 
airfoil presented herein will be based on the critical location at the leading 
edge adjacent to the platform as Indicated In figure 1, unless otherwise Indi- 
cated. This location was established from the finite element analysis by 
determining the Gaussian point which exhibited the largest total strain change 
during a mission cycle. There was some practical difficulty In determining 
this location because of the large number of elements, Gaussian Integration 
points and load-time Increments Involved and the consequent need to survey a 
vast amount of computer output printout. 

A number of cases were analyzed In which the material yield strength was 
progressively and arbitrarily reduced to Increase the severity of the plastic 
strain until reversed plastlcty was Induced. Comparisons are made between the 
stress-strain cycles computed from the simplified and MARC finite-element 
analyses. The comparisons are limited to the first mission cycle because of 
the exorbitant computing time required for the nonlinear finite-element analy- 
sis. As mentioned previously, creep analyses were not performed because the 
dwell times were too short for the temperatures Involved to have significant 
creep strains and the creep properties of the material were not adequately 
defined. 

Since the simplified procedure Is basically uniaxial, It can not directly 
account for material anisotropy. The most convenient assumption to have made 
was that the anisotropic effects could be neglected and the stress-strain his- 
tory approximated by using only the longitudinal properties. However, the 
questlonableness of this assumption Is Indicated In figure 5 which shows the 
difference In the computed stress-strain cycles between finite-element analyses 
using anisotropic and only longitudinal material properties. The results 
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presented In figure 5 were from MARC elastic-plastic analyses of the HPFTB 
airfoil In which the material yield strength was deliberately reduced for 
analytical purposes to Increase the severity of the cycle. The stress-strain 
hysteresis loops In figure 5 were for a different Gaussian Integration point 
than was subsequently determined as being the critical location based on the 
maximum cyclical total strain range criterion. 

The MARC stress-strain cycle calculated from the orthotropic properties 
of the directionally solidified MAR-M 246 alloy Is shown In figure 6. All of 
the plastic strain occurred during heatup on the preignition part of the mis- 
sion. The calculated plastic strain was small (under 300 microstrain) and 
confined to a local region at the leading edge. Initially, MARC elastic analy- 
ses were conducted for all 124 cycle load-time Increments used In the elastlc- 
piastlc analysis. Effective elastic moduli were obtained throughout the cycle 
from the computed effective elastic stresses and total strains. The simplified 
procedure was modified to use these effective elastic moduli to simulate the 
effects of the material anisotropy. Using this approximation, very good agree- 
ment was obtained between the MARC and simplified analytical cycles as shown In 
figure 6(a) . 

The problem was then rerun using only six elastic finite-element solutions 
as Input for the simplified analysis. These solutions were for the start and 
end points of the mission cycle and for the minimum and maximum temperature 
points during preignition and Ignition. To establish a more complete history 
of total strain, another 114 load-time Increments were obtained by Interpola- 
tion from the Initial elastic solutions. The computed stress-strain cycle 
from the simplified analysis using the reduced number of elastic solutions 
also shows reasonably good agreement with the MARC cycle (fig. 6(b)), although 
not quite as good as when the larger number of elastic finite-element solutions 
was used as Input. A noticeable discrepancy Is seen In the compressive strain 
region where the reduction In elastic analysis points resulted In failure to 
capture some of the cycle fluctuation due to transient thermal effects during 
the rapid engine cooldown. The CPU time for the six elastic finite-element 
analyses amounted to 5 percent of that required for one cycle of the nonlinear 
finite-element analysis. 

To Increase the severity of the problem, a series of analytical cases 
were run In which the material yield strength was arbitrarily and progressively 
reduced until the occurrence of plastic strain reversal on the unloading part 
of the cycle. The maximum plastic strain for this case was over 5000 micro- 
strain. Calculated stress-strain cycles from the two analytical methods are 
compared In figure 7. The same history of total strain as was created pre- 
viously from the 6 elastic finite-element solutions was used as Input for the 
simplified analysis for this case. The simplified analysis cycle In figure 7 
showed reasonably good agreement with the MARC results. Again the exception 
-was-a— 1 ow-stres s-r-egl on- durl ng uni oadl ng-where-the -severe-thermal— f-1 uc-t-uat-1 ons 
due to cooldown were not fully taken Into account with the reduced number of 
elastic analyses. However, the region of the cycle where this discrepancy 
occurred was elastic and would not have a significant effect on life prediction 
based on the calculated stress-strain response. The stress-strain cycle pre- 
dicted from the simplified method provided the stress/straln ranges and mean 
stress values normally needed for life prediction purposes to almost the same 
degree of accuracy as the nonlinear finite-element analysis. The CPU time per 
cycle for the simplified analysis was less than 0.01 percent of that required 
for the MARC elastic-plastic analysis. 
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SUMMARY OF RESULTS 


A simplified Inelastic procedure for calculating the local stress-strain 
history In a thermomechanlcally cycled structure was further developed to 
handle material anisotropy. This was accomplished by the use of effective 
elastic moduli that were determined from anisotropic finite-element analyses 
for a number of points In the mission cycle. The simplified analysis was 
exercised on airfoil problems for the first-stage high-pressure fuel turbopump 
blade of the space shuttle main engine. Predicted stress-strain cycles for 
the critical airfoil location were compared to stress-strain cycles computed 
from elastic-plastic finite-element analyses using the MARC code. The follow- 
ing general conclusions were drawn from the evaluation of the Improved simpli- 
fied procedure: 

1. The stress-strain response predicted from the simplified analysis was 
generally In very good agreement with the elastic-plastic finite-element solu- 
tions. The predicted stress-strain cycles provided the basic Information nor- 
mally needed for life prediction, such as stress and strain ranges and mean 
stress, to almost the same degree of accuracy as the finite-element analysis. 

2. Limiting the elastic finite-element analyses to several key points In 
the mission cycle and Interpolating between these solutions to create a more 
complete history of total strain, resulted In some Inaccuracy In Intermediate 
parts of the cycle due to the neglect of transient thermal fluctuations during 
the engine cooldown phase. However, the region where this discrepancy occurred 
was elastic and would not significantly affect the accuracy of life predictions 
based on the calculated local stress-strain response. 

3. The simplified procedure computed the stress-strain history at the 
critical location of the structure using about 0.01 percent of the CPU time 
required for MARC elastic-plastic finite-element analyses. There was an over- 
head computing cost for conducting elastic finite-element analyses of key 
points In the mission to define the Input total strain history. This addi- 
tional cost amounted to about 5 percent of the CPU time used In just one cycle 
of the MARC analyses. 
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TABLE I. - OS HAR-H 246 PHYSICAL PROPERTIES 
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TABLE II. - OS MAR-M 246 STRESS-STRAIN 


PROPERTIES (LONGITUDINAL) 
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Stress, 
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MPa 



21 °C 

649 °C 

816 °C 
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808 

875 

.2 

830 

855 

930 

.4 
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965 
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945 
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1 .0 

870 

960 

980 
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.00149 

.00156 
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(c) Blade rotational speed. 

Figure 2. - Mission cycle used for analysis. 
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Figure 3. - Airfoil finite element 
model. 
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Figure 4. - Airfoil temperature cycle. 
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Figure 5. - Comparison of Marc stress-strain cycles for critical location using 
anisotropic and longitudinal material properties. 
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124 ELASTIC FINITE -ELEMENT COMPILATIONS 



TOTAL MICROSTRAIN 

(a) 124 Elastic finite-element computations. 

Figure 6. - Comparison of simplified and Marc stress-strain cycles at critical location 
using DS Mar-M246 properties. 


6 ELASTIC FINITE-ELEMENT COMPUTATIONS 



(b) 6 Elastic finite-element computations. 
Figure 6. - Concluded. 
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